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ABSTRACT 

In this paper, we present explicit inverses for two Brownian-type matrices, 
which are defined as Hadamard products of certain aheady known matrices. The 
matrices under consideration are defined by 3n — 1 parameters and their lower- 
Hessenb erg-form inverses are expressed analytically in terms of these parameters. 
Such matrices are useful in the theory of digital signal processing and in the theory 
and applications of test matrices, i.e., matrices with known explicit inverses, 
which are thus appropriate for testing matrix inversion algorithms. 

Subject headings: Brownian matrix — explicit inverse — Hadamard product — 
Hessenberg matrix — numerical complexity — test matrix 



Introduction 



Herb old! (119691 ) gives the explicit inverse of a matrix G„ = {(3ij\ with elements 



i > j. 



Valvil (I1977I ) gives the explicit inverses of two symmetric matrices K 



and 



with elements 

respectively. K is a special case of Bro wnian-typ e matrix and On is a lower Brownian matrix 
as defined in ICover fc BarnettI Jl986h . Earlier, IPicinbonol Jl983h has used the term "pure 
Brown ian matrix" for the type of the matrix K; and lCarayannis. Kalouptsidis. fc Manolakis 
(119821 ) have treated the so-called "diagonal innovation matrices" (DIM), special cases of 
which are the matrices K and A^. 
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In the present paper, we consider two matrices Ai and A2 defined by 

Ai = K oGn and A2 = N o Gn, 
where the symbol o denotes the Hadamard product. These matrices have the exphcit forms 



and 





kih2 


kih 


klbn-1 


kibn 


kicii 




k2h 


k2bn-l 


k2bn 


kiQi 


k2a2 


ksbs 


ksbn-i 


ksbn 


kiQi 


k2a2 


k^az 


kfi-lbn-l 


kn-lbn 


kitti 


k2a2 


kzaz 


kn—lO"a—l 


h h 


kibi 


A;2&2 


ksbs 


kn—lbn—1 


knbn 


k2ai 




ksbs 


kn-lbn-1 


k'ab'n 


k^ai 


k^a2 


ksbs 


kn-lbn-l 


knbn 


'n-lC^l 


kn-lO'2 


kn-i^a 


kn—lbn-1 


k h 
'^n'^n 


kntti 


kna2 


kntts 


■ ■ ■ kfiQifi—i 


k h 



(1) 



(2) 



In the following sections, we give the explicit inverses and determinants of these matrices. 

2. The Inverse and Determinant of Ai 
2.1. The Inverse 

The inverse of Ai is a lower Hessenberg matrix expressed analytically by the 3n — 1 
parameters defining Ai. In particular, the inverse A^^ — [aij] of Ai has elements given by 
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the relations 



Ci-lCi ' 

kici' 

bn-l 



i-1 



aij = < 



n 

u=j-l 

I 



where 



I 0, 



i = j = 1, 



J - « > 1, 



i = 1,2, . . . ,n - 1, Co = 1, c„ = 6„. 



di — A;j_|_iaj+i6j — /Cjaj6j+i, i — 1, 2, . . . , n — 2, (io — ^i. 



— ^i) 
fl'i — ~ ki, 



i — 2,3, . . . ,n — 1, 



with 



i-l 



H kj^^i if 

and with the obvious assumptions 



(3) 



(4) 



(5) 



/ci 7^ and ci^Q, i — l,2,...,n. 



(6) 



2.2. The Proof 

In order to prove that the relations (3) give the inverse matrix A^^, we shall reduce 
Ai to the identity matrix by applying elementary row operations. Then the product of the 
corresponding elementary matrices gives the inverse matrix. In particular, we apply the 
following sequence of row operations. 



Operation 1. row i — x row (i + 1), i = 1, 2, . . . , n — 1, 

which gives the lower triangular matrix 

I {k2bi - hai) ... 

'^(h-k2) ^^(hb2-k2a2) ... 

^{k4-ks) ^{k,-ks) f^{k,h-k,a,) ... 



\^ {kn kn—l) {kn kn—l) (^n ^"n-l) ■ ■ ■ (^n^n— 1 ^n-lC^n-l) 

/ciOi A;2a2 /csOs ■ ■ ■ k 



Operation 2. row i — j^.'^fg. ^ x row (i — 1), i — n,n — 1, . . . ,3, kn+i — 1, 
which results in a bidiagonal matrix with main diagonal 

'kiCi k2C2 kn-lCn-l 



and lower first diagonal 

/ kiaig2 k2k^g^f2 kn-2kn-ign-lfn-2 kn^i 

V ka ' k4g2 ' ' A;„5(n-2 ' gn-i 



Operation 3. row 2 — ^jpM. x row 1, and 



row i — /cj — ^ ^.^^ _ i = 3, 4, . . . , n, 

i^i-\-l9i—l^i—l 



which gives the diagonal matrix 

kiCi k2C2 

k2 k-i 



kn—lCfi—l 
kn 



knCfi 
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Operation 4. x row i, i — 1,2, . . . ,n, 
which gives the identity matrix. 

Operations 1-4 transform the identity matrix to the following forms, respectively. 

Form 1. The upper bidiagonal matrix consisting of the main diagonal 

(1 , 1 , ... , 1) 



and the upper first diagonal 



ki k2 
k2 ' k3 



kn-l 



Form 2. The tridiagonal matrix 



1 








1 

fe4fl2 








^492 








1 + 







kn-2gn-l 
fcnSn-2 
fen 









1 + 



9n-l . 



Form 3. The lower Hessenberg matrix 



1 

_ 0132^2 

Q'lg3fc3fc2/2 
/C4C0C1C2 



fcl 

k2 



ksd 



fe4ClC2 



„ aign-lfc«-lfc2/2.--fcw-2/w-2 „ C^ign-lfcw-lfc3/3---fcn-2/n-2 
fenC0Cl...Cn-2 fc„ClC2...Cn-2 



^aignknk2f2---kn-lfn-l 



^dignknkzfz...kn-lfn-l 



COC\...C„-\ 



ClC2...Cn-l 



where the symbol s stands for the quantity (—1 



kji- 







{knbn-2 — kn-2(ln-2) 







kfi- 



h r o V'm'^n— ^ I'^n—A'-^n—Al u 

^n'^n — 2 n-n 
dn—2kn brt—lkn 



Cn-2Cn-l 



Cn-1 



Form 4. The matrix whose the elements are given by the expressions (3). 
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2.3. The Determinant 

The determinant of Ai takes the form 

det (Ai) = kibn (A;2&i - kiai) {kzh2 - A;2a2) ■ ■ ■ {kj)n-i - kn-ian-i) . (7) 
Evidently, Ai is singular iiki — or, using the relations (4), Ci — for some i e {l,2,...,n}. 



3. The Inverse and Determinant of An 



3.1. The Inverse 

In the case of A2, its inverse ^ = [cty] is a lower Hessenberg matrix with elements 
given by the relations 

kj-ibj-i — fcj-i-iaj-i 



where 



with 



a 



Ci-lCi ' 

J_ 

Cl' 

k-nCn—lCn ' 

i-1 



i = j = 1, 



-ly+L 



v=j+l 



n 

v=j-\ 



I 



I 0, 



di — kiai^ibi — /cj+iajftj+i, 
fi — (^i bi, 

9i — ~ 



J - « > 1, 



i = 1,2, . . . ,n - 1, Co 
i = 1,2, ...,71-2, do 
i = 2, 3, . . . , n — 1, 
i = 2,3, 1, 



1, 
ai. 



n 



1 if i = J + l, 



" =b 



(8) 



(9) 



(10) 



- 7- 



and with the obvious assumptions 

7^ and Q 7^ 0, i = 1, 2, . . . , n. 



3.2. The Proof 

In order to prove that the relations (8) give the inverse matrix ^2 ^) shall reduce 
A2 to the identity matrix by applying elementary row operations. Then the product of the 
corresponding elementary matrices gives the inverse matrix. In the present case, we carry 
out the following row operations. 



Operation 1. row i — row (i + 1), i = 1,2, . . . ,n — 1, 



which gives the lower triangular matrix 



kibi — k2ai 





A:2&2 - ^302 



ai {kn-i - kn) a2 (A;„-i - kn) 

kndl kn(l2 






kn-lbn-1 kndn-l 
kn^n—l 








h h 



Operation 2. row n — x row {n — 1) and 

row i — X row (i — 1), i = n — l,n — 2, . . . ,3, 

Qi-l V ' ) ) ) 

which results in a bidiagonal matrix with main diagonal 

(ci , C2 , ■ ■ ■ , c„_i , /c„c„) 

and lower first diagonal 

93k2f2 9n—lkn—2fn—2 ^n^n— l/n— 1 



am , , ■ ■ ■ , 

92 



9n-2 



9n-l 



Operation 3. row 2 — x row 1, 



row I 



0-192 
ci 

9iki-lfi-l 



X row (i - 1), i = 3, 4, . . . ,n - 1, and 
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row n - X row (n - 1), 

which results in the diagonal matrix 

C2 ... Cfi—i kfiCfi^ . 



Operation 4. ^ x row i, i = 1,2, . . . ,n — 1, and 
X row n 



Ci 

1 



which gives the identity matrix. 

Operations 1-4 transform the identity matrix to the following forms, respectively. 

Form 1. The bidiagonal matrix with main diagonal 

(1 , 1 , ... , 1 , 1) 

and upper first diagonal 

(-1 , -1 , ... , -1 , -1). 
Form 2. The tridiagonal matrix consisting of the main diagonal 
1,1, 1 + ^, ... , 1 + ^, 1+ ^" 



92 gn-2 Qn-l 

the upper first diagonal 

(-1 , -1 , ... , -1 , -1) 

and the lower first diagonal 

5*3 Qn-i k„ 



0, 



92 gn-2 9n-l 



Form 3. The lower Hessenberg matrix 

1 -1 ... 



_ aiS2 fci&i-fcaai 

CQCl Cl 

aig3fc2/2 dig3 

C0C1C2 C1C2 








„ nign-lfc2/2...fcn-2/n-2 „ C^lgn- 1 ^3 /s ■ ■ /n-2 fcn-2fcn-2 — fcraart-2 

C0C1...C„_2 Cl...C„_2 ■ ■ ■ Cn-2 

„ aignknk2f2---kn-lfn-l „ dignfonfes/s-.-fcrt-l/n-l dn-2gnkn kn-lbn-1 

CoClC2...C„_l ClC2...C„_l ■ ■ ■ Cn-2Cn-l C„_l 
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where the symbol s stands for (—1)*+-'. 

Form 4. The inverse matrix given by the formulae (6). 



3.3. The Determinant 

The determinant of A2 has the form 

det{A2) = Kbn {hbi - /caOi) (A;2&2 - ha2) ■ ■ ■ {K-ibn-i - Kan-i) , (12) 

which shows in turn that the matrix A2 is singular if A;„ = or, adopting the conventions 
(6), Cj = for some i G {1, 2, . . . , n}. 



4. Numerical Complexity 

The relations (3) and (8) lead to recurrence formulae, by which the inverses 
and A^^, respectively, are computed in O(n^) multiplications/divisions and 0{n) addi- 
tions/substractions. In fact, the recursive algorithm 



a. 



i,j+i = -l/ci, i = l,2,...,n- 1, 



bi-iQi -no ^2 bn-1 

aii^-ai^i+i-\ , z = 2, 3, . . . ,n - 1, an 



Ci—lCi fClCi Cn—lCn 

di-29i ■ r, o 

ttM-i = ' z = 2,3, 

Ci_2Ci_iCi 

_ dj-s-^kj-sfi-s . ri. -. r\ ■_9 

<^i,i— s— 1 — 7 (^i,i—sj ^ — O, 4, . . . , 71, S — 1, Z, . . . , Z Z, 

Oj— s— lCj_s_2 

where the q, dj, /j, and gi are given by the relations (4), computes A^^ in 5n^/2 + 5n/2 — 6 
mult/div (since the coefficients of ai^is depends only on the second subscript) and 5n — 9 
add/ sub. 

In terms of j, the above algorithm takes the following form 

ttj-ij = -1/ Cj-i, i = 2, 3, . . . , n, 

<^jj = > J = 2, 3, . . . , n - 1, an = 7 — , = , 

Cj—lCj KiCi Cji—lCji 
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dj-iQi+i -10 1 
= ' J = l,2,...,n- 1, 

a,+s+i,j = - ^^•+-+^^^+-^^+- aj+sj, j = l,2,...,n-2, s ^ 1,2, . . . ,n - j - 1. 

9j+sCj+s+l 

For the computation of ^2^^ the above algorithms change only in the estimation of the 
diagonal elements, for which we have 

an ^ -ai^i+i -\ , z = 2,3, 1, an ^ -a^, a^n = t , 

Ci—\Ci KjiC-n—lC-n 

where the q, rfj, /j, and (7^ are given by the relations (9). Therefore, considering the relations 
(4) and (9), it is clear that the number of mult/div and add/sub in computing A^^ is the 
same with that of A^^ . 
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